How the hydrogen bond in NH 4 F is revealed with Compton scattering 
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In order to probe electron wave functions involved in the bonding of NH4F, we have performed 
Compton scattering experiments in an oriented single crystal and in a powder. Ab initio calculations 
of the Compton profiles for NH4F and NH4CI are used to enlighten the nature of the bonds in the 
NH4F crystal. As a consequence, we are able to show significant charge transfer in the ammonium 
ion which is not observable using other methods. Our study provides a compelling proof for hydrogen 
bond formation in NH4F. 

PACS numbers: 78.70.Ck, 71.15.-m, 61.05.cf 



The distinction between ionic and hydrogen bonding 
is not always clear [![. Ammonium salts, in general are 
similar to the potassium and rubidium salts in crystal 
form and in other physical properties Q and thus are 
clearly ionic in nature. This similarity is due to the size 
of the ammonium ion, NH^, which is almost equal to 
the size of K + and Rb + ions. Remarkably, NH4F, in- 
stead of having a CsCl or NaCl structure like other am- 
monium salts which are ionic in nature, is isostructural 
to ice Ih with NHj and F~ ions alternating in the oxy- 
gen positions. Besides, ammonium fluoride and ice can 
form mixed crystals [3|. The fact that ammonium fluo- 
ride deviates from the standard isomorphism and prefers 
to adopt the ice structure is ascribed to the presence of 
significant hydrogen bonding in this compound 0, 0] ■ 

While ice has been the subject of many studies, NH 4 F 
has not been as widely studied. The crystal structure has 
been determined by Adrian and Feil [6( , the electrostatic 
properties have been studied by van Beek et al. 0, Ab 
initio calculations have been carried out by Alavi et al. 
[|| and by van Reeuwijk et al. [9]. Given the ionic nature 
of ammonium salts, it is tempting to consider the force 
between the NH|" and F~ ions to be essentially an elec- 
trostatic interaction. Nevertheless, charge density distri- 
bution does not provide a good criterion for establishing 
the existence of hydrogen bonds, as shown for example 
by Alavi et al. [8| for this compound. In this work we 
investigate the nature of chemical bonding in NH4F and 
show, using the study of electron momentum distribu- 
tions, that we detect the presence of hydrogen bonding 
in ammonium fluoride. 

Recently, Compton scattering, which is inelastic X- 
ray scattering at large energy and momentum transfers 
[lfj . [llT | has provided fundamental information on the 
quantum nature of the hy drog en bond in ice and wa- 
ter 0, El Q 0, [3, EE [Of- Thc measured Comp- 



(EMD) p(p): 



J(Pz) = j j P(P)dpxdpy, 



(1) 



ton profile (CP), J(p z ) represents the double integral of 
the ground-state total electron momentum distribution 



where p z lies along the scattering vector of the x-rays. 
Our Compton scattering experiments on NH4F will show 
how electronic states are modified by the hydrogen bond 
in contast to the ionic scenario of NH4CI. 

Thc CPs of the crystalline a, c directions and the spher- 
ical average of a powder were measured at the Euro- 
pean Synchrotron Radiation Facility in Grenoble on high 
energy beamline (ID15B), using the scanning Compton 
spectrometer [lij . The energy of incident photons was 
set to 29.74 keV by use of a Si(lll) cylindrically bent 
crystal and the scattering angle was 173°. The collected 
Compton spectra were energy analyzed using the Si (400) 
analyzer and detected by a Nal scintillator counter. Thc 
energy-dependent resolution function is deduced from thc 
full width at half maximum (FWHM) at the thermal dif- 
fuse scattering peak which was 0.16 (a.u.) which provides 
at the Compton peak 0.13 a.u. of resolution. After sub- 
tracting the background, energy dependent corrections 
have been performed: absorption in analyzer and detec- 
tor, detector efficiency, analyzer reflectivity [2(J. After 
normalization to the number of electrons, the contribu- 
tion due to multiple scattering was subtracted from the 
total measured profile [2l[ . In order to obtain the valence 
electrons CPs, a quasi-self-consistent-field calculated core 
CP is then subtracted from the total corrected measured 
profile 

The program used for our ab initio calculations was 
CRYSTAL98 [13], which is especially suitable for molec- 
ular crystals. We have used the hexagonal structure of 
the NH 4 F crystal at room temperature with point group 
P63mc, a = 4.37 A c = 7.17 A and the experimen- 
tal atomic positions [7|. This crystal structure is illus- 
trated in Fig. [TJ The occupied orbitals used to determine 
the EMD and the CPs were calculated using a restricted 
Hartree-Fock scheme [la] . 
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FIG. 1: (Color on-line) 48 atom cluster extracted from the 
NH4F crystal. The (red) largest spheres are F atoms, the 
small (blue) spheres are H atoms and the medium (green) 
spheres are N atoms. The crystallographic directions a, b and 
c are shown along the edges of a box containing the cluster. 



FIG. 3: Comparison of the theoretical spherical averaged va- 
lence CPs and the experiment for the powder. The dots are 
the experimental points. The dashed line is for the superpo- 
sition of neutral NH3 and HF molecules. The solid line is for 
the NH4F crystal. The curves are normalized to unit area. 
Symbol size is representative of error bars. 
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FIG. 2: Comparison of the theoretical and experimental [100] 
valence CPs. The dots are the experimental points. The 
dashed line is for the superposition of isolated NHj and F - 
ions. The solid line is for the NH4F crystal. The curves are 
normalized to unit area. Symbol size is representative of error 
bars. 



In Fig. [2 we show two computed valence CPs for NH4F 
(convoluted with experimental resolution and normalized 
to unit area) together with the experimental data along 
the crystallographic direction a. In the first calculation, 
represented by a dashed line, we have built the CP as 
a superposition of CPs of isolated ions NHj and F~ ar- 



ranged in the same geometry as in the crystal. In com- 
parison with the experiment, the profile of the isolated 
ions is found to be higher due to the absence of bonding 
17| . In the second calculation, represented by a solid 
line, the two ionic fragments can interact, and an excel- 
lent agreement between the crystal calculation and the 
experiment is obtained. 

Likewise, in Fig.[3l we show that the spherical averaged 
CP for the NH4F crystal calculation (solid line) agrees 
remarkably well with the experimental CP of the pow- 
der in contrast to a calculation for two non-interacting 
neutral molecules NH3 and HF (dashed line). In this 
molecular computation, the intramolecular atomic dis- 
tances are the same as the corresponding distances in the 
crystal. Interestingly, the spherical averaged CP for two 
non-interacting neutral molecules NH3 and HF turns out 
to be very similar to the spherical averaged CP for the 
isolated ions NH|" and F~. In reality, the main difference 
between the ions and the molecules cases is an exchange 
of a proton since the HF molecule can be formed by the 
capture of a proton by the F~ from the ammonium ion. 
However, the X-ray Compton scattering has not direct 
access to the proton wavefunction. 

In order to study bonding effects, we now discuss the 
CP anisotropy [l2j, that is, the difference between mea- 
sured Compton profiles along two crystalline directions. 
In the present case, we consider 



A{p z ) = Jc{Pz) - Ja{Pz), 



(2) 



where J c (p) is the CP along the crystallographic direc- 
tion c which is also the direction along which hydrogen 
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FIG. 4: Comparison of the theoretical and experimental 
Compton anisotropy [001] - [100]. The circles are the exper- 
imental points while the solid line is for the ab initio simula- 
tion of the NH4F crystal. The amplitude of the experimental 
anisotropy has been rescaled as in Ref. [13 ]. 



FIG. 5: (Color online) Pd(z) or the power spectrum of the 
Compton anisotropy for NH4F (experiment and theory) and 
NH4CI (theory). The curves are normalized to unit area. 



bonds arc aligned, and J a (p) is the CP along the crys- 
tallographic direction a. The anisotropy A(p z ) has the 
advantage of highlighting the changes introduced by the 
bonding. The shape of the anisotropy is very similar to 
that found in ice A(p z ) [12| as shown in Fig. 2] 

The effects produced by the hydrogen bond on the elec- 
tronic wavefunctions can be extracted by studying the 
power density 12. |25| 



Pd{z) 



dp z A{p z ) exp(ip z z) 



(3) 



The peaks in the power density Pd{z) indicate charac- 
teristic distances over which wave functions are coherent 
in given crystallographic directions. Figure [5] shows that 
the experimental data and the theory are in good agree- 
ment concerning the locations of the two main peaks. 
Nevertheless, the long range peak at about 4 A in the 
experimental anisotropy is largely reduced in the calcu- 
lation. A cause of this discrepancy could be the Hartree- 
Fock tendency to overestimate wavefunction localization. 
However, the CP in momentum space is not too sensitive 
to this overestimation. 

In a true ionic crystal the main differences from the 
isotropic picture occur in the region where the ionic 
spheres are in contact as in the hard sphere model. Thus, 
in the picture developed by Pattison et al. [26[ one ex- 
pects peaks at distance corresponding to the ionic diam- 
eters. The first main peak of Pd{z) should correspond to 
the cation diameter while the second main peak should 
give the anion diameter. Figure [5] also compares NH4F 
to the more ionic salt NH4CI [27J. One can see that the 
second main peaks for NH4F and NH4CI j2?| are close to 



accepted values for F - and Cl~ diameters (2.7 A and 3.6 
A respectively, 28]). The first peak is however at a dis- 
tance systematically smaller than that of the hard sphere 
ionic diameter of ammonium (2.8 A) even in NH4CI sig- 
nifying that the purely ionic, hard sphere model is not 
adequate. Moreover, the distance corresponding to this 
peak is smaller for both the theoretical calculation and 
the experiment in NH4F than the calculation for NH4CI 
though the ion is identical. As for NH4F, the interpre- 
tation becomes straightforward when the connection is 
made with results obtained for the hydrogen bonded ice 
Ih since the measured and calculated anisotropy bear a 
striking resemblance. In the more appropriate hydrogen 
bond picture the two predominant peaks of Pd{z) located 
at about 1.7 A and 2.75 A can be assigned to the hydro- 
gen bond length F-H (1.67 A) and the nearest neighbor 
F-N distance (2.77 A), respectively implying the exis- 
tence of coherent charge transfer as in the case of ice Ih 

In conclusion, our study shows that high-resolution di- 
rectional CPs of ammonium fluoride provide clear sig- 
natures of the hydrogen bond. The ionic picture is just 
a starting point since highly charged states such as the 
NH4" and F~ are mitigated by the hydrogen bond: bond- 
ing electrons are shared in states that cannot be assigned 
solely to ammonium or F. 
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